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We investigate the conditions under which a gap vanishes in the spectrum of dimerized coupled 
spin-1/2 chains by means of Abelian bosonization and Lanczos diagonalization techniques. Although 
both interchain (J') and dimerization (8) couplings favor a gapful phase, it is shown that a suitable 
choice of these interactions yields massless spin excitations. We also discuss the influence of different 
arrays of relative dimerization on the appearance of non-trivial magnetization plateaus. 
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Current developments in the field of spin ladders have 
revealed intriguing features of low-dimensional quantum 
antiferromagnets ( AF) [jl| . The appearance of plateaus 
in the magnetization curves of these systems has received 
much attention from both the theoretical and experimen- 
tal side ||]. A necessary quantization condition for the 
appearance of such plateaus in generic ladders was de- 
rived in Q (see also || for the case of ID chains and fl-fl] 
for other particular cases). In a very recent paper [pLan 
unexpected phenomenon has been observed in the mag- 
netization curves of NH4CUCI3 at high magnetic fields. 
The crystal structure of this material at room tempera- 
tures is known to be composed of double chains (2-leg lad- 
ders), with three different nearest neighbor interactions. 
Thus, from previous theoretical studies of such arrange 
of couplings in ladders systems, plateaus at (M) = and 
1/2 are expected to be visible at sufficiently low temper- 
ature. However, measurements performed in B at very 
low temperatures, clearly show two net plateaus at 1/4 
and 3/4, while no plateaus are observed at the theoreti- 
cally expected values (M) = and 1/2. 

As a first step to reconcile some of these facts with the 
current understanding of ladder materials, in this work 
we study two important issues namely, (i) the possibility 
of closing the gap in a two-leg dimerized ladder by a com- 
bined effect of the dimerization and the interchain cou- 
pling and, (ii) the emergence of a plateau at (M) = 1/2 
depending on the realization of the relative dimerization 
(see Figs. l(a)-(b)). Though (i) was studied on a quali- 
tative level using non- linear sigma model techniques , 
here we present a quantitative and more systematic treat- 
ment, whereas (ii) is analyzed for the first time in this 
work. These two issues could shed light in the study of 
the experimental measurements such as those performed 
in j|] as well as in related aspects of CuGeO chain com- 
pounds 1 10 1, in which dimerization becomes staggered be- 
tween weakly coupled chains. 

As is well known, the isotropic spin-1/2 Heisenberg 
chain is already in a critical state. Thus any relevant 
per turbation, such as the Peierls dimerization instability 
jjllj or a weak interchain coupling [jljl, [|13|| , can dras- 
tically alter the nature of the ground state, whereas a 
massive spin gap excitation appears simultaneously in 
the energy spectrum. Interestingly, it was suggested that 



a combined effect could lead to a massless regime ||. 

In the present work we study this issue quantitatively 
by means of both bosonization and numerical techniques, 
and find that there exists indeed a fine-tuning of the cou- 
plings responsible for the appearance of massless lines. 
Though this is a fine-tuning effect, our result implies 
that there is a finite region of the coupling parameter 
space where the gap is expected to be small. Thus, when 
the latter becomes comparable with thermal fluctuations, 
measurement attempts of zero magnetization plateaus 
could be smeared out even at low temperatures. 

In analyzing the abelian bosonization of spin-1/2 
Heisenberg ladders, we follow a similar methodology de- 
veloped as in |,|, Q, 0, Q, which is particu- 
larly suitable to elucidate the behavior of weak coupling 
regimes. Specifically, here we study two dimerized spin 
chains interacting through a Hamiltonian 
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(a = 1,2), where the S n denote spin-1/2 operators. For 
staggered ladders, e.g. Fig. 1(a), the array of coupling 

n+1 , and parametrized by 



exchanges are set as J„ = J„ 



Jn^ = J [1 + (— l) n <5] say for chain (1), whereas for the 
non-staggered situation we just set J„ = Jn , a s shown 
in Fig. 1(b). To maintain pure AF and non- frustrated 
exchanges, the dimerization parameter is kept bounded 
by 1 8 1 < 1 throughout the 2L spins of the ladder length 
with periodic boundary conditions (L even). 

On general grounds 0, it is expected that gapful mag- 
netic excitations should appear for all magnetizations 

(M) = i J2 n (Sn^ lS> + Sn ), satisfying the quantization 
condition pN(l — (M})/2 € Z, where p stands for the 
polymerization Q and for the number of chains cou- 
pled together. This would imply the existence of plateaus 
at (M) = and, if dimerization occurs (p = 2), also at 
(M) = 1/2, for a N = 2-leg ladder material such as the 
one studied in ||. However, the experimental magneti- 
zation curve shows no plateaus at these two values of the 
magnetization, but instead at (M) = 1/4 and (M) = 3/4 
||. To explain this apparent breakdown, we are espe- 
cially interested to include an homogeneous field term of 



the type - f £„ 



q z(2) 



so as to unravel the 
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interplay between the above kinds of coupling arrays and 
applied magnetic fields h, say along the z-direction. 

It is well known that the low-energy properties of the 
Heisenberg chain, (6 — 0), are described by a c = 1 con- 
formal field theory of a free bosonic field compactified at 
radius R = R({M),A) for any given mag- 

netization (M) and XXZ anisotropy |A| < 1 (see e.g. 
|jl6f ) . The functional dependence of R can be obtained 
using the exact Bethe Ansatz solution by solving a set 
of integral equations obtained in JT^Jl^l and [[l9) (for a 
fuller review consult for instance Ref. pj). The bosonized 
expression of the low-energy effective Hamiltonian for a 
single homogeneous chain, (8 = 0), in the presence of an 
external magnetic field h and with an XXZ anisotropy 
— 1 < A < 1, can be shown to adopt the form 

H=[dx~ {u 2 (x)/(4R) 2 + R 2 (d x cf>(x)) 2 } , (2) 

with II = d x <t>, and 4> = 4>l +4>r, 4> = 4>l~ 4>r- Here, the 
effect of the magnetic field h and the XXZ anisotropy 
A enter through the radius of compactification R. This 
radius governs the conformal dimensions, in particular 
the conformal dimension of a vertex operator e 1 "^ is given 

In the limit of both weak dimerization \S\ <C 1, and 
interchain coupling J'/ J -C 1, the bosonized action reads 

X 

A 2 cos{4k F x + V4^(0 (1) + <t> {2) )) + 

X 

A 3 ^ cos (^4^ (1) -<t> {2) )) + 

X 

A 4 ^cos(v^(0 (1) -0 (2) )) + 

x 

2 

^a^(-l) x cos{2k F (x + l/2) + ^^ a \x)). 

(a) = l X 

where Aj cx J'/ J, cx 6 and cx ±5, (the minus 
sign corresponding to Fig. la), and the Fermi momentum 
k F is related to the total magnetization (M) via kp = 
(1 — (M))n/2. In (j|) a marginal term has been neglected 
for the sake of simplicity since it does not change our 
results. 

(i) Closing of the gap 

At zero magnetization (i.e. kp = tt/2), all terms in (|^) 
are commensurate and relevant. Due to the Ai term, we 
have to first diagonalize the derivative part of the Hamil- 
tonian which is achieved by introducing new variables 
± = (cf) 1 ± <j) 2 )/\/2. The associated compactification 
radii are given by R± = R^Jl ± J 1 /(AJ^R 2 ). 

For AF interchain coupling, (J 7 > 0), and in the new 
basis, the A4 term is the most relevant and orders the 4>~ 
field (the dimerization terms mix the </> + and 4>~~ variables 



and will be treated as a perturbation). The a interac- 
tion in (j3[) is then wiped out after integrating the massive 
4>~ field, since it always contains a contribution from cj)~ . 
Thus we are left with the </> + field with a relevant inter- 
action given by A2. Hence, this field is in general massive 
so we can expect a plateau at zero magnetization. The 
interesting point here is that the a interaction generates 
radiatively a term that can cancel the A2 perturbation, 
(which is the one responsible for the mass of the sym- 
metric boson). From this plain weak coupling analysis 
we see that this happens on the critical line given by 

j*S 2 . (4) 

Therefore, whenever this cancellation occurs we have one 
massless degree of freedom corresponding to the symmet- 
ric combination of the original free bosons 4> + ■ 

To enable an independent check of this result, we now 
turn to a numerical finite-size analysis of the original lad- 
der Hamiltonian. We direct the reader's attention to Fig. 
2 where we display the energy gap necessary to create an 
excitation with total spin S = 1 in the staggered coupling 
array [Fig. 1(a)]. The results were obtained from exact 
diagonalizations of finite ladder lengths, 4 < L < 12, via 
a recursion type Lanczos algorithm poj . In studying the 
mass gap extrapolations towards their thermodynamic 
limits (dashed curve of Fig. 2), we fitted the whole set of 
finite-size results using a variety of standard procedures; 
namely linear, logarithmic and van den Broeck-Schwartz 
type methodologies of convergence Mj. Basically, they 
yield analogous results with at least 2 significant digits, 
though within the neighborhood of the massless region 
finite size effects become noticeable. In particular, this 
latter variation is reflected in the localization of the crit- 
ical line contained within the error bars of Fig. 3. Nev- 
ertheless, at least in the weak coupling (S, J') regime we 
were able to obtain a reasonable agreement with Eq. (Q). 
The progressive departure from the latter however, yields 
(3)he right critical parameters. As \S\ — > 1, we clearly re- 
cover the "snake chain" limit where criticality is achieved 
when J' = 2 J [see Fig. 1(a)]. This strong coupling regime 
yields in turn the expected linear behavior in (|<5| — 1) . 

Preliminary estimations of the sound velocity P2j near 
the ground state are consistent with a conformal anomaly 
c = 1 throughout the critical line, except for S = J' = 0, 
where the system has two decoupled bosons with c = 2. 
Moreover, a low temperature computation of the mag- 
netic susceptibility revealed that the latter is enhanced 
significantly on approaching the massless line (ji|), thus 
lending further support to our bosonization approach. 

Also, this effect is presumed to occur on the ferromag- 
netic side, J' < 0, where a critical line can emerge with 
usual dimerization as well [Fig. 1(b)]. In fact, a similar 
analysis using the perturbative RG approach has been 
given in the first reference of |l3| and, indeed, the effect 
is observed numerically [ ^3][ . 

(ii) Plateau at (M) = 1/2 
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Let us now consider Eq. (|3|) at non-zero magnetiza- 
tion (Af). As was referred to above, this can be read- 
ily accounted through the radius of compactification 
R((M), A) resulting from the effect of the magnetic field 
h and the Fermi momentum Uf = (1 — (M))ir/2, Then, 
the A2 and a interactions are incommensurate, and we 
have in principle a massless mode corresponding to the 
field. However an interesting phenomenon can be ob- 
served: a plateau in the magnetization curve is expected 
to appear at 1/2 of the saturation value due to the ap- 
pearance of the radiative correction 

L 

7^(-l) x cos(4Mx + l/2) + V8^> + ). (5) 

x=l 

This term is generated from a combined effect 
of the dimerization on each chain (a^ interac- 
tions in (|3|)) and the interchain interaction terms 
(Af) (cos(2k F x + \/47r</> (1) ) + cos(2fc F a; + V^(/> {2) )) ■ 

However, the presence of this radiatively generated 
term is strongly dependent of the manner in which the 
dimerization is realized. In the case of normal dimeriza- 
tion, one finds that 7 oc SJ'(M), while in the staggered 
case it vanishes exactly, 7 = 0, due to the relative minus 
sign between ceW and a^ 2 \ This effect, predicted from 
our weak coupling analysis can be clearly observed by 
comparing Figs. 4(a) and 4(b) (see below). 

The operator (^) is relevant for any anisotropy A > 0, 
where a Kosterlitz-Thouless (KT) transition is expected 
to take place. (More precisely, this KT point will depend 
also on the interchain coupling J' through R± 

One can also understand this difference in a sort of 
strong coupling limit, when \S\ — > 1. It can be read- 
ily checked that the normal dimerization displays an 
(Af) = 1/2 plateau Vh £ (2 J, 4 J) with J' = J, while 
there is no plateau in the staggered array. This effect is 
therefore robust in the sense that it is valid for all val- 
ues of the couplings and could then be a useful tool for 
experimentalists to decide whether the former or the 
latter exchange topology is actually realized in a given 
compound. 

To test the correctness of this scenario, finally we re- 
cur to the Lanczos algorithm, to obtain the lowest en- 
ergy eigenvalues on each magnetization subspace S z — 
{0, 1, ... ,L}. Using ladders up to L = 12, in Figs. 4(a) 
and (b) we display respectively the ground state mag- 
netization curves of the staggered and normal coupling 
arrays. By setting J'/J = 2 and <5 = 0.5, our numerical 
analysis supports entirely the bosonization picture as no 
plateau shows up for staggered arrays at any finite mag- 
netization [Fig. 4(a)]. In contrast, the normal array [Fig. 
4(b)] clearly exhibits the expected (Af) = 1/2 plateau. 
Furthermore, the robustness of this result was checked 
numerically for a variety of dimerization (5) and inter- 
chain (J') coupling regimes. In all cases, the staggered 
(normal) array always displays a massless (massive) be- 
havior at (Af) = 1/2. Finally, it is worth remarking that 



this study of the simplest TV = 2 leg ladder and with 
dimerization, (p = 2 in the notation of 0), can be ex- 
tended using the same bosonization methods to the study 
of N leg ladders with arbitrary polymerization p 0] . 

In summary, by using bosonization and numerical tech- 
niques, we have studied two interesting phenomena oc- 
curring in a dimerized two-leg ladder: (i) the closing of 
the gap at zero magnetization by means of a fine-tuning 
mechanism and, (ii) the disappearance of the 1/2 plateau 
(gap at finite (Af)) by alternating the dimerization of 
the chains along the rung direction. As for the first is- 
sue, we presented a thorough quantitative study of this 
phenomenon for AF systems, whereas (ii) was studied 
for the first time here. On this latter respect, though 
NH4CUCI3 under high magnetic fields is supposed to be 
described by a two-leg ladder system, its magnetization 
curves does not show up neither (Af) = nor (Af) = 1/2 
plateaus || . We trust the present paper will help to con- 
vey a better understanding of these ladder materials. 
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FIG. 2. Gap spectrum for different sizes of the staggered 
spin ladder for J'/J = 1. The dashed lines display extrapo- 
lations to the thermodynamic limit. 
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FIG. 1. Schematic view of alternating ladders with inter- 
chain coupling J' and dimerization parameter 8 for (a) stag- 
gered and, (b) non-staggered arrays 
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FIG. 3. Critical line of the staggered spin ladder in the 
(S, J') coupling parameter space, extrapolated from finite 
samples. Solid lines are guide to the eye whereas the dashed 
parabola stands for the bosonization approach. 
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FIG. 4. Magnetization curves 
8 = 0.5 and J'/J = 2 
(b) non-staggered coupling arrays. Solid, dashed and short 
dashed lines denote respectively magnetizations for L — 12,10 
and 8. The thick full line denotes the expected form in the 
thermodynamic limit 
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